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Abstract 

We examine the Ising chain in a transverse field at zero temperature from the point 
of view of a family of moment formalisms based upon the cumulant generating function, 
where we find exact solutions for the generating functions and cumulants at arbitrary 
couplings and hence for both the ordered and disordered phases of the model. In a t- 
expansion analysis, the exact Horn-Weinstein function E(t) has cuts along an infinite set 
of curves in the complex Jt-plane which are confined to the left-hand half-plane 3 Jt < 
— 1/4 for the phase containing the trial state (disordered), but are not so for the other 
phase (ordered). For finite couplings the expansion has a finite radius of convergence. 
Asymptotic forms for this function exhibit a crossover at the critical point, giving the 
excited state gap in the ground state sector for the disordered phase, and the first excited 
state gap in the ordered phase. Convergence of the ^-expansion with respect to truncation 
order is found in the disordered phase right up to the critical point, for both the ground 
state energy and the excited state gap. However convergence is found in only one of 
the Connected Moments Expansions (CMX), the CMX-LT, and the ground state energy 
shows convergence right to the critical point again, although to a limited accuracy. 
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1 Motivation 



The purpose of this work is to study the predictions of a family of non-perturbative moment 
formalisms devised to treat Field Theoretic and Many-Body Models with a well known exactly 
solvable model, the 1 — D Ising Model in a transverse Field. This model has a quantum phase 
transition, and we compare those predictions of the formalisms with the exact results in the 
two phases and near the critical point of the Model. These formalisms, the t-Expansion[[?J and 
the Connected Moments Expansion (CMX)||], are related in that they are centred upon the 
cumulant generating function, which takes an exponential form. These two Expansions have 
been applied in a number of areas ranging from quantum chemistry, lattice gauge theories 
to condensed matter problems and the reader is referred to Witte []T5| for references. The 
formalisms take their starting point from a finite sequence of low order connected Hamiltonian 
moments, or cumulants \v n : n = 1 . . . }, calculated with respect to a chosen trial state, with 
the proviso that the overlap with exact ground state is nonzero. In the i-Expansion, one 
constructs the cumulant generating function from its Taylor series expansion about t = 0, 

<e^>=exp{f>^}^), (1) 

where all moments are calculated with respect to this trial state. Then the Horn-Weinstein 
Theorem Q states that the ground state energy and first excited state gap (within the same 
sector as the ground state) || are given by the following limits 

E = lim E(t) 

t— >oo 

F F — lim d2E W dt2 (2) 
Ei-Eo-^So- dE{t)/dt > 

where E(t) = F'(—t) is termed the Horn-Weinstein function. Other ground state averages can 
be found by computing the appropriate generalised connected moments or higher excited 



states [11] found from certain derivatives of the Horn-Weinstein function and performing 
the same extrapolation. This theorem was "proved" in Ref. (R) by assuming a simple 
discrete point spectrum for the system, even in the thermodynamic limit, and these kinds 
of assumptions continue to be commonly used without adequate justification. The question 



immediately arises of how this extrapolation can be carried out and several methods [14] are 
employed. These extrapolation methods include Pade Analysis, where E(t) is assumed to be 
representable by a [L/M] Pade approximant with L = M i.e. it has only separable poles and 
zeros of finite order in any finite region of the complex t-plane, and the extrapolation t —* oo 
made directly, Other methods are D-Pade Analysis, Laplace Method/Resolvent Analysis, 
Inversion Analysis, and U E of F" /Partition Function Method, and the reader is referred to 
Witte|i~5|] for the original references. 

The Connected Moments Expansion (CMX) is founded upon the assumption of a par- 
ticular form for the asymptotic behaviour of E{t) as real t — > oo, 

n 

E^-Eo + Ys^e-^ , (3) 

j=i 

although its variants employ additional and differing assumptions. The variants give the 
following explicit forms for the ground state energy, according to their treatments. The first 
method, CMX-HW0], gives the ground state energy up the nth order as 



Sf 



E = v 1 -J2 ™ —~ > ( 4 ) 



i n}=i Sj,3 



2 



where the S-sequence are recursively generated by 



Si+ik 



SihSi 



s' 



Another method, CMX-SDQ, gives the ground state energy up the nth order as 

Eq = U>i >n , 

where the w-variables are recursively generated by 



(5) 



(6) 



^j+l,k^2,k 



^3,k 



Uj+2,k 



(7) 



The last method, the CMX-LT described in Ref. 



gives 



Eq = V\ - (f2,^3 



/ ^3 ^4 
l/ 4 U 5 



"3 



, ■ ■ ■ , "nj 



V n +2 



(8) 



V Vn j 



called the AMX[10], which is very similar to the CMX-HW 



There is also another variant, 
method. 

Very little has been known of an exact nature concerning these methods. For example 
it was not known how the formalisms could be realised in an exact manner, whether the 
exact theorems were in fact confirmed for solvable Models, whether the assumptions made in 
the inevitable extrapolations needed to treat non-integrable Models were justified, and what 
factors affected the convergence of the methods. Recently such questions were begun to be 
addressed with respect to another solvable Model, but that Model was studied only at its 



critical point - the isotropic XY Model in l-Dimension[15|. We extend this analysis here to 
a closely related Model, within the same universality class, and reexamine these questions 
with regard to a phase transition. We firstly analytically solve our Model from the moment 
point of view, then investigate the validity of the theorems using the exact Horn-Weinstein 
function, following this we reveal the analytic character of this function in the complex t- 
plane, and in the truncated expressions using a set of finite low order cumulants examine the 
convergence of ground state spectral quantities with respect to the order of truncation. In 



a companion work to this paper, paper II[16], we continue analysis of this model from the 
point of view of the other family of moment methods - the geometrical family, which includes 
the Lanczos Method. 
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2 The ITF Model 



The model we investigate is the Ising model in a transverse field in one dimension, which 
exhibits a order-disorder quantum phase transition at zero temperature^, M. The Hamil- 
tonian for this Model is taken to be the following 



N N 
n=l n=l 

which differs slightly in conventions with previous use. We define a coupling constant ratio 
x = B/J. In what follows the label n indexes the sites and periodic boundary conditions 
fj^r +1 = af are applied. This model has a zero-temperature quantum phase transition at 
B = J or x = 1, which separates an ordered phase x < 1 whereby the z-component of all 
the spins are aligned either up or down from the disordered phase x > 1 induced by the 
tunnelling term. This model is also self-dual, in that H(x) — ► xH(l/x) leads to the same 
model, and consequently for the ground state properties. The exact ground state is given by 



eo /J = -h 1 + x \E( J-^-) , (10) 

where E{k) is the complete Elliptic function of the second kind, according to the conventions 
of Gradshteyn and Ryzhik||]. The first excited state gap has the exact form 

AE = 2J\l-x\ , (11) 

which vanishes at the transition point. For the full details of the solution to this model and 
discussion of its properties we refer the reader to the works [12, |J. 



In the first step to an explicit form for the cumulant generating function one makes a 
transformation from the ordered spin operators to the disordered spin operators 

n-l 

= Y\. a m ' = a n a n-l ) = MnA*n > (12) 

m=l 

and then the Wigner-Jordan transformation to a pure Fermionic representation 

In doing so we have neglected the boundary terms arising from the exact mapping of the spin 
or a-cyclic problem, and will only consider the c-cyclic problem cn+i = c\. This followed by 
a Fourier transformation of the Fermionic operators 



N 

1 



^I>^c n , (14) 

n=l 

where the momenta are k = 2nl/N and the integers I = —N/2, . . . ,N/2 — 1 for N even, 
whereas I = —(N — l)/2, . . . , (N — l)/2 for odd. Proceeding in this way the Hamiltonian 
becomes partially decoupled as 

H = h k = ^2 {2iJ sin k[cl c± k +c k c^ k ] + 2(J cosfc+J3)[cjj7%+cl fc c_ fc -l]} . (15) 

fc>0 fc>0 

For each k there is a 4-dimensional Hilbert space 1 00) , |01), |10), |11) but because c£ c^. , c k c_ k 
only connect 1 00) , |11) and ctck,c_ k c—k are diagonal one only needs the 2-dimensional BCS 
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reduced subspace |00),|11). Thus we have a pseudo-spin problem | JJ.) = |00),| ft) = |11) 
with pseudo-spin operators 

T k = 4 c ~k > T k = c -kCk , r| = 4c fc +c+ fc c_fe-l , (16) 

and thus h k = — 2Jsinfcr| + 2(Jcos k+B)r%. 

As an obvious first choice for the trial state, but not the only choice, we take the exact 
eigenstate as x — ► oo, that is the ferromagnetically aligned state \J-) polarised along the 
x-axis | JJ-) = | III • • • |||)z, which is represented in the fermionic c— picture by 1 000 . . . 000). 
Another equally obvious choice for the trial state would be the exact eigenstate for x = 0, 
and this we denote and is the fully polarised state in the z-direction. This implies that 

we only need the matrix elements with respect to 

(U \e th *\ U) = cosh(tE k ) - 2(JCO " fc+jB) smh(tE k ) , (17) 

where the quasi-particle energy is 

El = 4 J 2 + AB 2 + 8 JB cos k . (18) 
The moment generating function is then quite simply given by 

(e tH ) = II {cosh(t^) - 2{JC ° Sk+B) smHtE k )\ 

fJV f 1 ,,, [ . . 2(Jcosk+B) W „J1 
exp < — / dk In cosn.(tE k ) smh(tE k ) } ■ 

I 27r JO L E k J J 



iV-K» 



(19) 



To simplify our discussion we move to dimensionless quantities and the dimensionless quasi- 
particle energy squared E k = 2Je k is then 

e\ = l + x 2 + 2xcos/c . (20) 

It should be noted that it was not necessary to go any further and fully diagonalise the 
problem in order to find the cumulant generating function. 

However continuing with the diagonalisation process one employs a Bogoliubov trans- 
formation to quasi-particle operators r) q ,r)£ 



Cq = U q r] q - iv q r]± q , C- q = U q f]- q + W q r]+ 

4 = u i4 + il >qV-q , c± q = u q rj± q - iv q rj q 



(21) 



with transformation elements 



. e q + x + cosq\ 1/2 f e q -x-cos q \ 1/2 



2e ? J " { 2e q 



The trial state can now be expressed in terms of quasi-particle operators and the exact ground 
state by 



-,-r \l e +x + cos q + i\l e a — x — cos q rfirr ' 
\F) = II — %= ^=i|0) , (23) 

q >o y Ze i 



and the overlap squared simply reduces to 



N f w f e q +x+cosq 



HW r__y ^.i^-^]. (24) 
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As x —>■ oo then |(jF|0)| 2 — » 1, as it must, while for x = 1 we have 



ln|(.F|0) | 2 = (-ln2 + 2G/vr)iV 



(25) 



where G is Catalan's constant, and finally at x = 0, |(.F|0)| 2 = 2"^. Clearly the logarithm 
of the overlap squared is simply proportional to N with a negative coefficient which is zero 
at x = co and decreases monotonically to — In 2 at x = 0. 

The lowest excited state can be found from the moment problem associated with a trial 
state with the correct quantum numbers, and the obvious candidate is 



which can be shown to be proportional to The relevant moment generating function 

is then simply 




n=l 



(26) 




(£\e- tH \£) 



e -2t\J-B\ 



(27) 




because the quasiparticles are non-interacting in this model. 
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3 The Analytic Character of E(t) 



The analytic character of the Horn-Weinstein function E(t) in the complex t-plane is of fun- 
damental importance in devising any extrapolation strategy for the t-expansion, as knowledge 
even of a qualitative nature or the lack of it is the greatest source of errors in extrapolation. 
From the moment generating function the exact Horn-Weinstein function E(t) is given by 

vr J 1+ x+cosfc tanh(2Jt£fc) 

Before treating the case of general x, the two extreme cases of x — > 0, oo are taken first, 
which lead to simple results. The limiting form of E(t) when x — > oo is, of course, completely 
trivial, namely E{t) — > —NB. As will become apparent later however the way this limit is 
approached exhibits some complex character. In the case when x — > some care in the limit 
needs to be taken and the following result holds 

m= {-NJt a nh(Jt) (4tf-l)f O(Jt) < (4tf+l)f ReZ (2Q) 

\-NJ coth(Jt) (4K"-3)f < 3f(Jt) < (4*T-l)f 

The analytic form this function has is one with an infinite number of Riemann sheets bounded 
by cuts at Jt = (2K+l)ji for all K 6 Z. There are no divergences in this function anywhere, 
while it has discontinuities on the branch lines defined above. Correspondingly the radius of 
convergence of the expansion about t = is not vr/2, as might be expected from considerations 
of the growth in the cumulants, but is it J A. 

The singularity structure of E(t) for general x is most interesting as it is a singular 
integral, and has discontinuities across a denumerable set of open, infinite curves in the 
complex t-plane. This set of singular curves, where the denominator of the integrand of 
Eq. (28) vanishes, is defined by 

3f Jt= (K+l/2)^- K€Z 

1 , ( 1 /1-x 2 M ( 3 °) 

$Ut= — tanh-M — 



2e {2x 

and parametrised by a real e with \x — 1| <e<x + l when 1, and by 

QJt = (K + 1/2)^ KeZ 
4£ 

1 , ( 31 ) 

Jt = tanh" 1 ^) , 

when x = 1, with < £ < 1. This family of curves in the complex Jt-plane can be simply 
described as the different Riemann sheets of the following transcendental equation 

\/x + Z — y/x + l/z 



exp [2JtVxTzy / x+ 1/z) = v v ; 1 , (32) 

v > ^x + z + Vx + 1/z 

with an arbitrary complex parameter \z\ = 1, and where the only choices of the branch or 
sheet arise from the logarithm and not the square-roots. In Figure (0) these curves are shown 
for a sequence of x-values from 4 to 0, and a remarkable change occurs at the transition point. 



7 



For x > 1 the singular curves are confined to the left-half of the complex Jt-plane, to the left 
of 



RJtr 



cosh 1 (x) 



1 



(33) 



with both asymptotes of the curves bent back parallel to the real axis. At x = 1 the 
asymptotes furthest from the real axis have opened up and become parallel to the imaginary 
axis. They are still confined to 3ft Jt < —1/4. However once x < 1 these asymptotes invade 
the right hand side of the Ji-plane, becoming parallel to the real axis as K Jt — ► +oo, although 
tending towards a different 9 Jt value from that of the left-hand asymptote. Finally at x = 
the singular curves are straight and parallel to the real axis. Thus we see that the big shift 
from one phase to the other is marked by the loss of analyticity in the right-hand side of 
the complex Jt-plane. Expressed in another way one can perform a Wick rotation from the 
positive real axis through ±7r/2 without encountering any singularities for x > 1, but for 
x < 1 this is not possible. With this information we can find the radius of convergence of the 
Taylor series expansion about t = 0, from the distance of the singular curves for K = 0, — 1 
from the origin. If we denote this by R(x), we have 



R(x) 



4e 



'^ 2 + In 2 



l-(e -xf 



(34) 



where the parameter eo is given the solution of 

l-(eo-x) 2 4xe (eo - x 2 + 1) 



In" 



(e +x) 



1 



+ 



[(x+l) 2 -e 2 ][e 2 -(x-l) 2 ] 



In 



1 - {e -xf 



(e +xy 







(35) 



The radius of convergence approaches tt /4 for small x-values and decreases monotonically 
as x is increased, finally taking the form of ln(2x)/2x for large x. This behaviour is shown 
in Figure (Q) . However the actual size of the radius of convergence is not important in this 
method, although it should be non-zero, but rather whether the right half of the complex 
Jt-plane is free of singularities. 

As a final step in describing the analytic nature of this function we give the following 
result for the discontinuity in it across the singular curves, 



AE(Jt) 



8iNJe 3 



Ve 2 - (x-1) 2 

2\ 



4xe(e 2 + l 



-x 



[(x+1) 



e 2 ]\e 2 



[x- 



1 



+ ln 



1 



+ m{2K + T 



(36) 



Asymptotically along either arm of the curves e— >|1 — x|,l + x this jump vanishes for all 
x values, while there are single, local but finite enhancements at some intermediate value 
of e (a slightly different point for the real and imaginary parts). This appears as a sharp 
enhancement for all curves which occurs close to the imaginary axis, but the jump here is not 



divergent (it is not possible for the denominator of Eq.36 to vanish). Thus all of the singular 
nature resides in discontinuity of the function in passing from one Riemann sheet to another. 
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4 The Asymptotic Forms of E(t) 



Of central importance in understanding the t-expansion theorems in our Model is the asymp- 
totic form for E(t) as t — > oo in a sector which contains the positive real axis. There several 
reasons for this. Firstly any changes in the asymptotic behaviour with coupling reflects 
changes in the analytic character of E(t) with coupling also, secondly it has a bearing on 
extrapolation strategies for the t-expansion, and lastly it gives the precise forms that should 
be assumed for Connected Moment Expansions. The explicit form for the Horn-Weinstein 
function satisfies the ground state Horn-Weinstein theorem in that lirm ; _ >00 E(t) — * eo where 
the ground state energy is given by Eq. ([H]). To find higher order asymptotic terms we have 
to consider the cases x > 1, x = 1, and x < 1 separately as we will find a crossover in the 
asymptotic behaviour at the transition point. In the first case, x > 1, one can make a simple 
expansion of the denominator of the integrand, integrate term-by-term, and apply a variant 
of a lemma of ErdelyiQ to arrive at 

k=l \ n=k L v /J 

1 + ' L> nmJtkx(x+l)]"+V2 nk[ } 



n=k 



NJ (x-l)V 2 4Jt(^i) 

7rV2 [8Jte] 3 /2 



(37) 



for |Arg Jt\ < 7r/2. This expression is not valid at x = 1. The U \ (x) are polynomials in x 
with integer coefficients, defined by the generating functions 



~( U nk(x) = ^{l + 2XZ) 2 [1 - X(x-I)z] k -V 2 [l + (x-l)z]- fc - 1/2 [l + XZY^I 2 , 

V Z (38) 

V ^UL(x) = z k (l + 2xz) 2 [l + x(x+i)z] k -V 2 [l + (x+l)z]- k -V 2 [l + xzY k - x l 2 . 



n! 

ra=0 



Clearly, in this phase, the asymptotic form has a dominant exponential decay controlled by 
the first excited state gap 4J(x — 1) in the ground state sector together with the algebraic 
factors t~ n ~ 1 l 2 for n > 1. There are also subdominant terms with an exponential decay 
controlled by the gap at the top of the spectrum AJ(x + 1), which are discontinuous as Arg Jt 
changes sign as a manifestation of the Stoke's phenomena. The Stoke's line is Arg Jt = 0, on 
which the first asymptotic series has maximal exponential dominance over the second series, 
or its associated series. The Anti-Stoke's line are Arg Jt = ±7r/2 where the asymptotic and 
its associated series are of the same order. For the purposes of the Horn-Weinstein theorem 
within this sector of the t-plane the dominant term in the asymptotic expression above is 
the exponential term with a unit multiple of the gap 4J(x — 1), then the higher algebraic 
terms combined with this, while all other exponential terms are subdominant. The relative 
ordering of these latter terms can vary with coupling but has no physical significance for 
the Horn-Weinstein theorems. It can also be shown that the excited state Horn-Weinstein 
theorem when applied to the above function E(t) yields 
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which arises from the excitation of two noninteracting quasi-particles near momenta k = 

7T + 0(l/iV). 

Exactly at the critical point the leading exponential factors vanish leaving the dominant 
series having only algebraic decay of the form, 



8NJ " 



E{t) - ~ -^l j2(-) k \ ^Jtkr^Vn 

r(n+l/2)r(l/2-Jfe)(n-2fc+l/2) ) 



'nk 

k=l Kn=l 



_l_ e -8Jtfc 



(n-k)\T(3/2-n)[16Jtk] n + 1 / 2 J 



with V^ifc defined by the generating function 

oo n 

£ Vn^) = *(1 - z) k -^(l + . (41) 

z — ' n! 

n=0 

This expression is still valid for |Arg Jt\ < vr/2, as before. The prediction of the Horn- 
Weinstein theorems from this case would be that the first excited state gap would vanish, 
and correctly so. 

However when x < 1 something interesting and different happens from that found for 
x > 1. Firstly it appears that an asymptotic expansion in the above sense is no longer 
possible, in that such an expansion should be valid within a finite sector containing the 
positive real axis. From the previous section we know that the Horn-Weinstein function is 
only analytic within a sector of zero extent about the positive real axis, and any asymptotic 
expansion can only be valid within this. Furthermore it is clear from the defining £-integrals, 



Eqj28, that for x < 1 the integrand has an inverse square-root singularity at the lower limit 
of the integral in contrast to the case of x > 1 which has a square-root branch point at 
both limits of integration. It is the algebraic singularities at the end-points of the integration 
t-interval which control the exponential factors in the asymptotic expansions. Instead of a 
full asymptotic expansion we recover just the leading order term, 

E{t) - ~ 2NJ(l-xfe- 2Jt ^ ) , (42) 

which implies that the excited state gap is now 

lim^^ = 2J(l-x). (43) 
t^oo E (t) 

In this phase the Horn-Weinstein excited state theorem gives the lowest excited state, rather 
than the next one lying in the ground state sector, because the trial state appropriate for 
x > 1 has a finite projection onto the first excited state for x < 1 (it is exactly orthogonal to 
the first excited state for x > 1, \£ )), and it gives the smaller of the two gaps. 

The lowest excited state gap is found from the generating function given in Eq. (|27| ) and 
it is trivially and exactly true that 

for all x which is just the energy of the lowest quasiparticle energy at k = n. 
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5 The Cumulants 



In the application of the t-expansion to non-integrable models, one computes cumulants with 
respect to the trial state up to some finite order and uses the truncated cumulant generating 
function in the ensuing analysis. Because our convergence studies will do just this, we need 
to compute the cumulants. Furthermore issues concerning the radius of convergence of the 
Taylor expansion of E(t) about the origin will arise in any discussion of the analytic character 
of this function. It is sometimes convenient to discuss normalised cumulants c n /J n which 
are polynomials in x and we will often use the same symbol to denote these as there can be 
no risk of confusion. The first few normalised cumulants are then 



ci = — x 

C2 = 1 

C3 = 4x 

c 4 = 16x 2 - 2 

c 5 = 64x 3 - 48x 

c 6 = 256x 4 - 672x 2 + 16 . 

The general form for the cumulant coefficients has an integral representation 
(2J) n+1 



(45) 



Cra+l 



dy 



2-7T y/ cosh 2 y — x 2 dy n 



tanhy 



y cosh 2 y — x 2 — x sinh y 
cosh y 



n+2 



found using the change of variables 



tanhy 



cos k+x 



(46) 



(47) 



The cumulant is a polynomial in x despite appearances which is a consequence of choosing 
the trial state to be an exact eigenstate of the second term in the Hamiltonian, Eq. (||). 
The explicit forms for these cumulants are then 



C2n+1 



(2J) 



2n+l 



(2p ! (n-p)l 1-x 2 n+1 .„ , , 
> — r 7— ——7-7- — -. — -J— — A™{ sgn(x 2 -l)(n+l+p)pP 



p=0 



l+x z 



1-x 2 



+ {n+l-p)Pf 



n+l 



1+x 2 



1-X 2 



C2n+2 



(2J) 



2rt+2 



^ (2p )! (n+l-p)!|l-x 2 r +1 p 
Zj, 2 P+ip! (n+l+p)! (2x)p p n+1 



1+x 2 



1-x 2 



(48) 



where the P™(z) are the associated Legendre Polynomials and the A™, are integer coeffi- 
cients defined from the higher derivatives of tanh(z), 



d 



2n 



dz 2n 

d 2n+1 
dz 2n + l 



tanh(z) = tanh(z) ^pSech 2p (z) , 



p=0 



n+1 



(49) 



p=0 
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One of the simplest explicit forms of these coefficients is one in terms of the Stirling 
numbers of the second kind, S(n, k), in the following way 

2n / 7 -i \ 

A n = { _ )P2 2(n-p) J2(-) k k\i^ p P _ 1 jS(2n,k), (50) 



v ' 



and the other coefficients are given by A™ + 1 = —2pB, 

The cumulants with respect to the trial state \J-&) are very simply given by 

c#(J,B) = c n (B,J) (51) 

and these normalised cumulants are then related by 

c*(x) = x n c n (l/x) , (52) 

i.e. the reciprocal polynomials. The cumulants simplify for the three special cases x = oo, 1,0. 
In the case x —* oo, where the J — > limit is taken and B fixed, then only the first cumulant 
is non-zero, c\ = —B, as would be expected. At the transition point x = 1 then 

c n+1 = (-4J) n+1 - dysechytanh n+1 y— tanhy , (53) 

which does not appear to be expressible in terms of any known mathematical constants, and 
forms a new integer sequence |lj|. At x = the odd cumulants vanish while the even ones 
are given by 

c 2k = (2J) 2k { 2k > B 2k . (54) 
where the B 2 k are the Bernoulli numbers. 
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6 Convergence Studies 



Of great practical interest is the question of convergence of the the various formalisms arising 
out of the i-expansion and how they change crossing the transition point, and here we examine 
all the dependences of the errors on the order of truncation for the ground state energy and 
the excited state gap as functions of coupling. In the results reported here all quantities were 
calculated using exact rational arithmetic including cumulants, truncated generating func- 
tions, Pade approximants, matrix inner products, solutions of the CMX recurrence relations, 
and CMX energies using the computer algebra system Maple V. Conversion to approximate 
floating point numbers (number of digits was 32) was only made for the final answer, so loss 
of precision or roundoff issues do not arise here. All data files and Maple source code are 
available upon request to one of the authors (NSW) . 

In the first study we employ only the simplest extrapolation strategy for the t-expansion, 
namely diagonal Pade approximation. In Figure @ we have plotted the relative error of the 
ground state energy evaluated by a [r/r] diagonal Pade approximation versus the truncation 
order r for various couplings x. Clearly the error is not smooth or strictly monotonic with r, 
as the Pade Approximant can change significantly from one order to the next but the overall 
trends are obvious. The general convergence is quite good, decreasing exponentially rapidly 
with order, with the rate of decrease highest for large x. As x approaches unity this rate of 
convergence drops, but doesn't vanish altogether, so we have the remarkable result that even 
at the critical point we still have substantial convergence. For couplings less than unity, such 
as the example given, we do not see any convergence beyond a minimal value. 

For the estimates of the excited state gap we expanded the derivative ratio, given in 
Eq. (|2|), as a Taylor series then formed the diagonal [r/r] Pade Approximations to this and 
made the extrapolation. In Figure (||) we have plotted this gap versus the truncation order 
r for the various couplings. Again for small orders there are wild fluctuations as the Pade 
Approximant shifts from one order to the next but these die out and converge quickly to 
the correct values, given by Eq (11). And again, most remarkably, the gaps converge as the 
coupling tends to unity, even at the transition point, and there there is evidence of a late 
convergence even in the other phase, at x = 1/2. 

In the CMX studies we first consider the CMX-LT variant. In Figure (|5|) we have plotted 
the relative error of the expression, Eq. (||) , truncated to r x r matrix versus the truncation 
order r for the same sequence of couplings x. While there are still fluctuations in going from 
order to order these are less than for the t-expansion, and the rate of convergence is much 
faster. The slowing down of the rate of convergence as the coupling approaches the critical 
one is also clearly evident, although even at the transition point there is still a respectable 
level of convergence. In the other phase, at x = 1/2 there is some improvement, even if 
the fluctuations are much larger, and the error saturates again at a modest value. For the 
HW variant of the CMX, the ground state energy estimate failed to converge to the correct 
value even at x = 2, rather it converged to another value with a relative error of 0.027%. 
In the case of the CMX-SD method it actually diverged for this coupling value, showing no 
convergent tendencies within the first 20 cumulants. 
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7 Conclusions 



In this work we have found the exact generating function and cumulants for the Ising Model 
in a transverse field for arbitrary coupling constant, with respect to a simple trial state 
taken to be the exact eigenstate at infinite coupling. For couplings larger than the critical 
coupling we have found that the Horn-Weinstein function is analytic in the right half of 
the complex plane, and has cuts along an infinite number of curves lying in the left half of 
the complex plane. When the coupling strength is less than the critical one, these singular 
curves cross both halves of the complex plane and the function is only analytic in a strip 
along the real axis. The radius of convergence of the Taylor expansion for the H-W function 
is finite, but does disappear for large coupling, however this is not quantitatively important 
for the t-expansion or Connected Moments Expansion. The Asymptotic form of the H-W 
function for large real positive argument displays a crossover as the coupling strength goes 
through the critical point. For couplings larger than the critical value one has a dominant 
exponential decay governed by the excited state gap within the ground state sector, along 
with an algebraic series of non-integral powers. At the critical point there is pure algebraic 
decay with integral exponents. And for couplings less than the critical one there is exponential 
decay determined by the first excited state gap. Results of the convergence for the t-expansion 
reveal that with the simplest extrapolation method, the ground state and excited state gap 
estimates converge in a systematic way for all couplings larger than the critical one, and 
even at the critical point, although with limited accuracy there. Only the CMX-LT variant 
converged, of the CMX methods, and in this case the ground state energies converged quite 
rapidly for all couplings in the above range, again even at the critical point. Because the 
actual character of the H-W function in these Models is rather different from that commonly 
assumed, and that these Models are not atypical of other strongly coupled Models, then 
extrapolation strategies in the t-expansion that employed this qualitative knowledge would 
be expected to have better convergence performance. Furthermore these examples provide 
concrete forms for new variants of the Connected Moment Expansion to be developed. 
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